Flat directions and gravitino production in SUSY models 
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Flat directions in supersymmetric models can get large vacuum expectation values in the early 
Universe which leads to a large mass for gauge bosons and gauginos. We point out that this can then 
result in enhanced gravitino production because the cross-section for the production of the ±1/2 
helicity states of the gravitino is proportional to the square of the gaugino masses. We consider 
gravitino production after inflation in such a scenario and find that the abundance in some cases can 
be much larger than the upper bound on the gravtino abundance from cosmological constraints unless 
the flat direction field has a very small vacuum expectation value when it commences oscillating. 

PACS numbers: 98.80.-k, 11.30.Pb 



I. FLAT DIRECTIONS IN SUSY MODELS 

Generic supersymmetric (SUSY) models have a large 
number of flat directions i.e. directions in the field space 
of scalars that have a flat potential [Tj-Q . However SUSY 
breaking and non-renormalizable terms lift these flat di- 
rections and can give rise to a non-zero minimum for 
the potential of the field <j) that parametrises a flat di- 
rection. During inflation the minimum of the potential 
for <j) can be large. <j> oscillates about its potential mini- 
mum, which varies with the Hubble parameter H . When 
H ~ mo, where too is of the order of SUSY breaking 
scale ~ 100 GeV , the minimum shifts to the origin with 
<f> displaced from the origin. In certain other cases, the 
minimum of the potential during inflation is at the ori- 
gin yet the flat direction vev is non-zero due to quantum 
fluctuations of <f> during inflation. 

The flat direction vevs give mass to gauge bosons 
through terms like g 2 A ll f*A fl f where / represents a 
squark or a slepton associated with <f>. The gauge bosons 
get a SUSY conserving mass m 2 ~ g 2 (p 2 where (p 2 = 
(<j>*<j>). In this process one or more Standard Model (SM) 
gauge symmetries can get broken (at late times <f> de- 
cays and gauge symmetry is restored.) Gauginos get an 



equivalent mass, and we take m g 



\foLLp. 



There can be two possible scenarios in the context of 
flat directions with large vevs. 

1. If all the gauge bosons get large masses, then ther- 
malization of the inflaton decay products is delayed. 
This is because gauge bosons mediate the scatter- 
ing processes relevant for thermalisation and the 
scattering cross-section is suppressed due to their 
large masses. Lack of thermalization then affects 
the cosmology of the Universe. In particular, one 
needs to recalculate the gravitino production for a 
non-thermal Universe, and for a lower final reheat 
temperature. This idea was first proposed in Rcfs. 



9] 



A more detailed calculation is done in Rcf. 



2. If not all gauge bosons get mass then thermaliza- 
tion happens as usual. However, the large vev of 
a flat direction gives a large mass to one or more 
gauginos which can still affect gravitino production. 



The gravitino has ±1/2 and ±3/2 helicity states and 
the total gravitino production cross section is propor- 
tional to ^1 + m~/(3m~)Y where the factor "1" is due 



to the ±3/2 helicity states and the factor to|/(3to~) is 
due to the ±1/2 helicity states [l(j. A large gaugino mass 
can lead to a large gravitino abundance in both the above 
scenarios. This has been overlooked in the literature so 
far. In this article we consider the second scenario in 
which only gluons and gluinos get mass and the Universe 
does thermalize quickly after inflation. The first scenario 
of a non-thermal Universe after inflation and enhanced 
gravitino production has been studied in Ref. Q. 



II. GRAVITINO PRODUCTION 

Gravitinos are produced by the scattering of the decay 
products of the inflaton. Refs. [U, [l2| provide a list of 
processes for gravitino production. The number density 
of gravitinos generated is then obtained using the Boltz- 
mann equation 



dt G 



(EtotM) n 2 



(1) 
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where n = (C(3)/7r 2 )T 3 is the number density of inflaton 
decay products (C(3) = 1.20206.. is the Ricmann zcta 
function of 3), E tot is the total scattering cross section 
for gravitino production, v is the relative velocity of the 
incoming particles, and (...) refers to thermal averaging. 
We ignore the gravitino decay term above as the gravitino 
lifetime is 10 7 ~ 8 (100 GeV /mg)s [ll[ and is not relevant 
during the gravitino production era for gravitinos of mass 
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10 2 3 GcV that we consider. (E to tM) is given by [0 



M 2 

1 3tt 

M 2 16 C(3) 



E 



3m| 



c i9i ln ( — 



(2) 



The i = 1,2, 3 in r](T) refer to the three gauge groups 
U(1) Y ,SU(2) L and SU(3) C respectively. g t {T) arc the 
gauge coupling constants. The parameters Ci and ki are 
constants associated with the gauge groups and are given 
by ci, 2 ,3 = 11,27,72 and ki, 2 ,3 = 1.266,1.312,1.271 re- 
spectively (see Table 1 of Ref. [3). M = Mpt/V&x ^ 
2.4 x 10 18 GeV is the reduced Planck mass. The cross- 
section above includes corrections to the expressions ob- 
tained earlier in Refs. [hJ EH- We presume that the in- 
flaton decays perturbativcly and the products thermalise 
quickly as discussed in Appendix A of Ref. [l5j . Wc 
shall only consider i = 3 because we are studying the 
enhancement in the gravitino production rate due to the 
large gluino mass and we take = 6 x 10 -2 for the 
relevant energy scales below. Then r](T) is given by 



3vr 



16C(3) 



3m l, 



(3) 



where rrig is now the gluino mass. 

To account for the expansion of the Universe one 
should consider the abundance of a species i in a comov- 
ing volume. This is achieved by considering the ratio 
Yi = tii/s, where rii is the number density of particles of 
the species i in a physical volume and s is the entropy 
density given by s = (27r 2 /45) g*T 3 , where = 228.75 
in the MSSM. Thus Eq. §j} is rewritten as 



. dY~ 

T^- = (X tot \v\)Yn. 



(4) 



In the radiation dominated era, the temperature T is 
given by 



T = T h 



1 



1/2 



2H R {t-t R ) + l 
where Tr is the reheat temperature and 



Br. 



8tt 3 5 



*R 



T 2 R 



90 Mpi 



(5) 



(6) 



This implies that T is 



Hr 



Thus 



dY, 



dT 



G _ 



90 



rpo 



1/2 



90 



1 



V2 T 3 
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(7) 



(2^ 2 /45) 5 , 



V_\ /C(3)V 

MJ \ 7T 2 J 



(8) 

To obtain the gravitino abundance we integrate the 
above equation from T R to Tf where / corresponds to 
the time when the flat direction condensate decays. Here 
we have ignored gravitino production during the period 
of reheating, as studied in Refs. [13, ITa - [20l |. Regard- 
ing tf, the perturbative condensate decay rate is g iven 
by T = ml /if 2 [J, [U and as argued in Ref. [2l[ the 
condensate decay products do not dominate the Universe 
for (p < 10 _2 Mp; for gravitational decay of the inflaton. 
(We discuss alternate mechanisms for condensate decay 
below.) For gravitino production after tf there will be no 
enhancement due to a large gluino mass and the abun- 
dance generated after tf will be proportional to Tf as 
similar to the standard scenario. 

We take the inflaton mass m^ to be 10 13 GeV and the 
inflaton decay rate Td ~ ?n 3 /Af|, ; ~ 10 GeV. Below wc 
are primarily concerned with the evolution of <p after to ~ 
m,Q . In the both the cases mentioned at the beginning 
of Sec. ((T| <fi effectively has a quadratic potential after to 
with a minimum at the origin with a positive curvature of 
TOq (ignoring thermal effects for now). The gluino mass 
is given by m| = aip 2 , and tp is tpo at to and then falls as 

1/a 3 / 2 once the condensate starts oscillating at to- Then 
for t > td, where td = T^ 1 is the inflaton decay time, the 
gluino mass is 
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mo 



T 
Tr~ 



(9) 



where we have used a ~ t 2 l z for to < t < td for an inflaton 
oscillating in a quadratic potential during reheating and 
a - l/T for t > t d . 



Integrating the Boltzmann equation in Eq. ([5]) from T R to Tf gives the gravitino abundance generated between td 



3 



and tf as 




(10) 



We have ignored the variation of and gi with temperature and have used the value at Tr, since for Tr 3> Tf most 
gravitino production occurs close to Tr. The term proportional to Tr — Tf is associated with the production of ±3/2 
helicity gravitinos and the other term is associated with the production of the ±1/2 helicity states. 



The reheat temperature is given by Tr = 

0. 55g* 1/4 (M Pl r d ) 1/2 = 2 x 10 9 GeV H3- The conden- 
sate decays when its decay rate T = m^/tp 2 equals H, 

1. e., when tf = (p 2 /rriQ, where 



a 

<Pf = <Po \ — 

■ a f 



3/2 



<P0 



3/2 / x 3/2 
CIO \ I &d 



a d 



3/4 



Here we have used a ~ t 2 / 3 for to < t < td and a 
for t > td- This implies 

2/5 1/2^1/10 

iff = ip Q m Q T d . 
The temperature at tf, using a ~ t 1 ! 2 for i > td, ii 
'ad' 



(11) 

, t l/2 

(12) 



= T R 



1/2 



2/5 -p3/5 



= 0.55 



-i/4 mpMpj 

2/5 t^I/10 



(13) 



For the condensate to decay after inflaton decay, or Tf < 
Tr, we need ipo > 30 mo or 93/ > 3 mo, for Yd = 10 GeV . 



III. RESULTS AND CONCLUSION 

A large gravitino abundance can be in conflict with 
cosmologies] observations. If the gravitino is stable its 
energy density can overdose the Universe. If the grav- 
itino is unstable its decay products can either overdose 
the Universe or dissociate the light nuclei generated dur- 
ing primordial nucleosynthesis. We use an upper bound 



of Yq < 10~ 12 for consistency with cosmological obser- 
vations. For the abundance obtained in Eq. QlOp to be 
consistent with observations then requires 



< 300. 



For mg = 100 GeV, the constraint is 



ip < 3 x 10 4 GeV 



(14) 



(15) 



Note that if we fix mg = 100 GeV and allow the con- 
densate mass, so far taken to be mo, to vary one gets 
the same constraint as in Eq. Q3] with mg replaced by 
the condensate mass (since in Eq. (1101) m|/ 

This is relevant if one includes thermal 
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corrections to the condensate potential. 

We obtain this bound assuming Tf <C Tr in Eq. ([TC 
From Eq. JT3J) tp = 3 x 10 4 GeV implies T f = 6 x 
10 8 GeV . For this value of Tf and T R = 2x 10 9 GeV , our 
assumption only slightly affects the contribution of the 
first term in Eq. (|T0|) . and does not alter the contribution 
of the second term. Since the contribution of the second 
term is larger than that of the first term our assumption is 
therefore justified. Smaller values of (fo will correspond 
to larger Tf and therefore greater compatibility of Yq 
with cosmological constraints. 

We now consider plausible values of ipo. During infla- 
tion the non-zero vacuum energy breaks SUSY and can 
give large positive masses to the flat direction of order 
Hi, where Hj is the Hubble parameter during inflation 
0- This may be true in supcrgravity models with min- 
imal or non- minimal Kahler potentials [23[ . The flat di- 
rection potential has a minimum at the origin and the 
field is driven to this minimum. Quantum fluctuations 
during inflation then give a vev of order Hi [H . Assuming 
that the field does not vary much till to, then ipo ~ Hi 
and the upper bound on ipo implies that the inflationary 
scale Vj /A < 4 x 10 n GcV. 

Alternatively, in no scale supergravity or more gener- 
ally in any supergravity theory with a Heisenberg sym- 
metry of the kinetic function in the Kahler potential, 
one gets contributions during inflation to the flat di- 
rection potential from supersymmetry breaking (due to 
Hi) only at the one loop level and these are negative 
at the origin for flat directions which do not involve a 
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stop and if the inflaton is not the standard string dilaton 
[23| . This correction to the potential, along with non- 
rcnomalizable terms, leads to a shifted minimum of the 
potential and subsequently a large vev of order Mpi [23| 
or 10 12 ~ 14 GeV on including GUT interactions [24[ as 
the field rolls to the minimum of its potential and oscil- 
lates about it. (Negative corrections at the origin were 
also discussed in Ref. 0].) 

We consider the shifted minimum of the potential for 
4> to be M( J ff/M) 1 / n+1 where non-renormalizablc terms 
in the potential are of the form cj) 2n+4 /M 2n , n > 1 (HH. 
4> oscillates about this time dependent minimum which 
decreases as H decreases. When H ~ mo, the potential 
minimum goes to zero and the field oscillates about the 
origin in a quadratic potential with curvature m§. Then 

a 

~ MiH^/M) 1 /^ , (16) 

where H(t ) = 100 GeV. For n = 1, (p < 3 X 10 4 GeV 
is satisfied for 

il/<10 7 GeV. (17) 

For larger values of n, the upper bound on M reduces 
further and one requires M < 10 4-5 GeV for n > 1. 
However we are considering scales upto at least Tr ~ 
10 9 GeV, and setting M = 10 9 GeV gives a value of tpo 
in conflict with the upper bound on ifQ. If we choose 
tfo ~ 10 12 GeV as mentioned above then we get the grav- 
itino abundance Y~ to be 8 x 10 2 which is orders of mag- 
nitude higher than the cosmological bound of 10~ 12 . The 
abundance will increase as ip\ for larger values of (fio. 

Thermal corrections to the condensate pot ential can 
increase the mass of the condensate field |25l - [27| . Ref. 
[Ill obtains masses of order 10 10 GeV, 10 6 GeV and 
10 5 GcV for non-renormalizable terms with n = 1,2,3 
and M = 10 18 GeV . We find that the condition 
(pa < 300 discussed above is not satisfied for ipo ob- 
tained from Eq. (|16[) with H(to) replaced by m^. 



Our results imply that either one has a scenario of 
supergravity with a positive contribution to the </> po- 
tential from Hj that gives Hi < 3 x 10 4 GeV and 
Vj^ 4 < 4 x 10 n GeV, or the flat directions must de- 
cay quickly in a way so as to not affect the cosmology of 
the Universe. 

With regards to the quick decay of the flat directions, 
the longevity of flat directions has been debated in Refs. 
[2ll. [28i- 36] . However it has been argued in Refs. (33l.l36| 
that even if non-perturbative rapid decay via parametric 
resonance occurs for scenarios with multiple flat direc- 
tions it leads to a redistribution of energy of the conden- 
sate amongst the fields in the D flat superspace and hence 
to practically the same cosmological consequences, in- 
cluding at least as large gauge boson and gaugino masses 
as in the scenario with only perturbative decay. 

Scattering of particles of the thermal bath off the flat 
direction condensate can lead to the decay of the conden- 
sate 0, H|| [2(| , though thermal effects are less important 
for larger values of n . For example, for n — 3 the con- 
densate decays much after the decay of the inflaton [2(| . 
Decay via fragmentation into solitonic states called Cu- 
bans [37i - [45l ] or Q-axitons [46| due to inhomogeneities in 
the condensate may also be relevant. However the time 
scales for the formation of Q-balls and Q-axitons can be 
larger than tj. tf is less than 60too > f° r <Pf < 800 GeV 
obtained from Eq. (|T2"j) . 

In conclusion, the presence of flat directions in SUSY 
models and associated large vevs seems to have very 
important consequences for gravitino production in the 
early Universe. These large vevs give a large mass to 
gauginos which enhances gravitino production and can 
violate cosmological constraints on the gravitino abun- 
dance. This scenario can be avoided if the flat direction 
has a vev smaller than 3 x 10 4 GeV at t ~ 1 , or if the 
flat direction decays early. 
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